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In the last decade, enormous progress has been achieved in the understanding of the 
various facets of coalescing double neutron star and neutron black hole binary systems. 
One hopes that the mergers of such compact binaries can be routinely detected with 
the advanced versions of the ground-based gravitational wave detector facilities, maybe 
as early as in 2016. From the theoretical side, there has also been mounting evidence 
that compact binary mergers could be major sources of heavy elements and these ideas 
have gained recent observational support from the detection of an event that has been 
interpreted as a “macronova”, an electromagnetic transient powered by freshly produced, 
radioactively decaying heavy elements. In addition, compact binaries are the most plau¬ 
sible triggers of short gamma-ray bursts (sGRBs) and the last decade has witnessed 
the first detection of a sGRB afterglow and subsequent observations have delivered a 
wealth of information on the environments in which such bursts occur. To date, com¬ 
pact binary mergers can naturally explain most —though not all— of the observed sGRB 
properties. This article reviews major recent developments in various areas related to 
compact binary mergers. 

Keywords: gravitational waves; nucleosynthesis; gamma-ray bursts. 
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1. Introduction 

In July 1974 Russel Hulse and Joseph Taylor discovered the pulsar PSR 1913-1-16 
during a systematic search for new pulsars at the Arecibo Observatory in Puerto 
Ricof^ The pulsar was detected with a period of 59 ms, but with apparent changes 
of 80 fis from day to day. Soon it became clear that these changes are due to the 
Doppler shift caused by the orbital motion around an unseen companion star. The 
orbital period is only 7.75 hours and with a semi-major axes of about 3 Rq, but 
no observed eclipses, the companion star must also be a compact object, either a 
neutron star (ns) of a black hole (bh). The binary system is so tight that the orbital 
velocity is of order 10“^c and general-relativistic effects build up qickly. To quantify 
the deviation of the dynamics from a purely Newtonian binary, one introduces so- 
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2 S. Rosswog 

called “Post-Keplerian” (PK) parameters. Since these PK parameters each depend 
on both component masses, one can determine individual stellar masses through 
the measurements of at least two PK parameters and not just -as in Newtonian 
theory- the total binary mass via of Kepler’s third law. PSR 1913+16 has by now 
been observed for more than four decades and its parameters have been determined 
to an astonishing precision. For example, the current valued for the neutron star 
masses are rup = 1.4398+0.0002 Mq for the pulsar and rric = 1.3886+0.0002 Mq for 
the companion star. The binary is a “clean” system in the sense that the components 
can be very accurately described as point masses and in particular tidal effects can 
be safely neglected. The periastron advances at a rat^of + = 4.226598 + 0.000005 
deg/yr (as compared to 43 arcsec per century for Mercury) and the orbit decays at a 
rate that agrees with General Relativity’s prediction for the emission of gravitational 
wave^ 

■ _ 1927rG^/3 1 + (73/24)e^ + (37/96)e^ 

^ 5c^ (TOp + TOc)^/^ V-Pfc/ {1 — e'^y/'^ 

to about 0.3%.!^ The latter, however, is only true if the relative Galactic accelera¬ 
tion difference between the pulsar and the solar system -they have a separatiorP of 
about (9.9+ 3.3) kpc- is taken into account. General Relativity had been probed in 
a number of solar system tests such as the perihelion shift of mercury,IS gravitational 
light deflectiorI30 or the Shapiro delay,or, more recently the Lense-Thirring effect 
in satellite orbitsjlEl but neither had the predicted radiative properties of gravity 
been probed nor regimes where the deviation from a flat space is substantially larger 
than in the solar system (gravitational potential $/c^ ^ 10“®). The orbital decay 
of PSR 1913+16 was the first -though still indirect- confirmation of the existence 
of gravitational waves and it laid to rest a longstanding debat£3IIi] about the exis¬ 
tence of gravitational waves and its quadrupolar nature. In 1993 Hulse and Taylor 
were awarded the Nobel Prize in Physics for their discovery of PSR 1913+16. 

In 2003 an even better suited laboratory for strong field gravitjElE3 was discov- 
eredEdEH the “double pulsar” PSR J0737-3039 where both components are ac¬ 
tive pulsars, usually referred to as pulsar A and B. With an orbital period of 
only 147 minutes, relativistic effects are even larger and with a distance of only 
about 1 kpc the systematic errors due to the relative Galactic acceleration are 
much smaller. The double pulsar has an eccentricity of e = 0.088 and shows with 
u! = 16.8991 + 0.0001 deg/yr an even larger periastron advance than PSR 1913+16. 
Its masses have again been very accurately measured as tua = (1.3381 + 0.0007) 
Mq and ms = (1.2489 + 0.0007) MQ^f^and within less than 10 years this binary has 
become the best testbed to probe the accuracy of gravitational quadrupole emission 
with a current accuracy well below 0.1%.!^^ 

To date over 2000 pulsars are known, about 10% of them possess a binary compan- 
ioiP^ and 10 systems consist of two neutron stars.l^^ Such systems have turned out 
to be precious laboratories to probe gravitational theories.li^t^ Although binary 
evolution suggests that it would be natural to form also a reasonable fraction of 
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neutron star black hole binaries pHUl to date no such system has been found. They 
should, in a number of respects, have properties similar to double neutron star bi¬ 
naries and we will therefore follow the common practice to collectively refer to both 
types as “compact binary systems”. 

Typically, one distinguishes gravity tests in different regimes: 1) quasi-stationary 
weak field, i.e. velocities v c and the spacetime is close to Minkowskian every¬ 
where, 2) quasi-stationary strong-field, i.e. velocities u <C c, but the spacetime can 
show substantial deviations from Minkowskian, 3) highly dynamical strong-field, i.e. 
V becoming comparable to c and involving a strongly curved spacetime and finally 
the 4) radiation regime where the radiative properties of gravity are tested. The first 
regime includes solar system tests, the second and fourth regime can be probed by 
well-separated neutron star binaries while regime 3) will be probed by the ground- 
based gravitational wave detectors such as LIGO, VIRGO and KAGRA.I^^EH They 
are expected to see the last inspiral stages minutes) that are initially well de¬ 
scribed by post-Newtonian methodd23211 and the subsequent merger and ringdown 
phase where strong-held gravity, hydrodynamic and nuclear equation of state (EOS) 
effects from the neutron star matter determine the dynamics. For an overview over 
the numerical modeling of these phases we refer to recent reviews 
Gompact binary mergers had also been suggested very early on as the production 
site for “rapid-neutron capture” or “r-process” elements,l2SHini actually before the 
discovery of the Hulse-Taylor pulsar. Although it may seem a natural idea to form 
r-process elements in the decompression of the extremely neutron-rich neutron star 
material, this idea has for a long time been considered as somewhat exotic. This 
perception, however, has changed during the last one and a half decades, mainly 
because compact binary mergers (CBM) were found to robustly produce r-process 
elements without hne-tuning while core-collapse supernovae seem seriously chal¬ 
lenged to providing the conditions for at least the heaviest r-process elements. This 
issue will be discussed in more detail in Sec. [H The idea has been further boosted 
by the recent discovery of an nIR transient in the aftermath of the short GRB 
130603B. Such transients from radioactive decays of freshly produced r-process ele¬ 
ments had been a prediction of the compact binary merger model. In particular the 
characteristic time scale of ^ one week and a peak in the nIR are consistent with 
extremely heavy elements with large opacities having been produced. This topic 
will be discussed in more detail in Sec. 14.41 

The merger of a compact binary system is also the most likely trigger of short 
gamma-ray burst (sGRBs). The idea that the sGRB phenomenon is caused by 
compact binary mergers has now stood the test of three decades of observation and 
it is considered the most plausible model for the “engine” behind sGRBs, although 
it is not completely free of tension with some observations. We will review in Sec. 
the major arguments that link sGRBs with the idea of an compact binary merger 
origin, we will discuss how the predictions from the compact binary merger model 
compare to the observed GRB properties and where open questions remain. 
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In the following, various facets of compact binary mergers are discussed in more 
detail. We begin by discussing relevant time scales in Sec.[^ followed by the expected 
neutrino emission in Sec.[^ nucleosynthesis and closely related issues are treated in 
Sec. 1^ and we discuss short GRBs in Sec. We will conclude with a summary in 
Sec. [6l 

2. Time scales 

As a stellar binary system revolves around its centre of mass it emits energy and 
angular momentum via gravitational waves, which are both extracted from the or¬ 
bital motior0with a power P oc ■ Since the total, kinetic plus potential, orbital 
energy is Aorb = —Gmim 2 / 2 a a loss of energy means a reduction of the mutual 
separation, which, due to Worb = \/GM/a^, corresponds to an increase of the or¬ 
bital frequency and therefore to a further enhanced gravitational wave emission. 
This runaway process finally leads to a merger of the two stars. For elliptical orbits, 
not only the semi-major axis a, but also the eccentricity e evolves in time. For a 
non-zero eccentricity, its the temporal change is negative, de/dt < 0, and the ec¬ 
centricity is “radiated away” very efhciently. The Hulse-Taylor pulsar, for example, 
has currently an eccentricity of e = 0.617 and a semi-major axis a = 2.2 x 10® cm. 
By the time the orbit has shrunk to 10 Rns, however, it will havd^ an eccentricity 
of e < 10“®, i.e. the orbit will be very close circular. 

A non-zero orbital eccentricity can substantially shorten the inspiral time until 
coalescence. In general, the evolution equations for a and e need to be integrated 
numerically, but a good approximation for not too large eccentricitiet|^ is given 

row 9,83 X lOVar* (]^) (j^) <2) 

Contours of the inspiral time log(TGw) for an equal mass binary system with 1.4 
Mq per star are shown in Fig. The inset shows the inspiral time (in seconds) for 
the last stages of a circular binary system as a function of the separation oq. Also 
marked are the locations of PSR 1913-1-16 and of the double pulsar PSR J0737- 
3039 although their masses are slightly different from 1.4 Mq. Circular binaries can 
only merge within the lifetime of galaxy (10^® yrs) if they possess an initial orbital 
period below 16 hrs. The Hulse-Taylor pulsar has a life time of 3 x 10® years until 
coalescence. A binary system with the same properties, but an initial eccentricity 
of e = 0.99, in contrast, would merge within only 10^ years. This dependence of 
the inspiral duration on eccentricity also has important implications for compact 
binary mergers as possible sources of r-process elements, see Sec. |4.3[ and for their 


®'Here, we will assume that we are dealing with a binary system of two non-spinning point masses 
in quadrupole approximation as first worked out by Peters and Mathews in 1963!^ 

'^This estimate ignores an additional eccentricity-dependent factor which is to within a few percent 
equal to unity for e < 0.6, and approaches monotonically 5^ 1.8 in the limit of e ^ 1, so that the 
estimates are always substantially better than a factor of 2. 
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Fig. 1. Inspiral time for a binary with 2 X 1.4 Mq. Shown are contours of log(TGw) (in years) 
as function of the current orbital period (in hours) and eccentricity. The filled circles indicate 
period and eccentricity of the Hulse-Taylor Pulsar PSR 1913+16 and the Double Pulsar J0737- 
039 (although their masses are not exactly 1.4 Mq). The inset focuses on the last inspiral stages 
where the orbits are to very high accuracy circular (note that the axes are different from the main 
plot). The last stages of the inspiral (uq < 100 km) only take fractions of a second. 


role as central engines of gamma-ray bursts, see Sec. 

The dynamics of the inspiral is also responsible for compact binary systems not 
being spun up to corotation when the coalescence occurs. For example, it takes only 
0.76 s for a 2 X 1.4 Mq binary system from the last 120 km of separation (« 10i?ns) 
until coalescence. Therefore, the stage where tides can be excited is simply too 
short to enforce corotation, even for an implausibly large internal viscosity 1^1^ 
Therefore, the spin period at merger is much smaller than the orbital frequency 
and an irrotational binary configuration is a sensible initial condition for modeling 
a binary system prior to merger. 

The frequency of the emitted gravitational waves /gw as a function of the time 
until coalescence tqw iP 


/1.21Mq 
V -^chirp 


5/8 


— ) 

TGW/ 


3/8 


/gw = 134Hz 


(3) 
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where the “chirp mass” is Mchhp = Therefore, a Mchirp = 1-21 Mq binary 

would emit gravitational waves of 10 Hz or more (roughly the lowest frequencies 
accessible to ground-based detectors) during the last 17 minutes and 100 Hz or more 
during the last two seconds of the inspiral. For a double neutron star system, one 
expects the final merger to occur around gravitational wave frequencies of 1 kHz. 
The inspiral can be accurately described by means of Post-Newtonian methods jUlllll 
but the final merger and “ringdown” require numerical simulations.^^^^^ Simulations 
indicate that the central object that forms in a merger does in many cases not 
collapse directly to a black hole, but instead survives as a “hypermassive neutron 
star” (HMNS),I2^1® despite having a mass in excess of the Tolman-Oppenheimer- 
Volkoff maximum mass. This is due to effects such as thermal pressure or -more 
importantly- differential rotatioi|^ A typical velocity field inside the merging binary 
system is shown in Fig.[^ at contact a shear interface forms for irrotational binaries 


^Differential rotation is much more efficient in stabilizing stars than uniform rotation since -in the 
center— the rotation can be very rapid without shedding mass at the surface. 




Fig. 2. Velocity fields during the merger of two initially non-spinning 1.4 Mq neutron stars.^^ 
The shear interface, see panel one, becomes Kelvin-Helmholtz unstable and forms a string of 
vortex rolls. Subsequently the vortex rolls merge and the remnant finally forms a radidly differ¬ 
entially rotating, so-called “hypermassive neutron star” that is temporarily stabilized against the 
gravitational collapse to a black hole. 
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Multi-messenger picture of compact binary mergers 7 

(left panel) that becomes Kelvin-Helmholtz unstable (right panel) and evolves into 
a differentially rotating remnant. Differential rotation is very efficient in stabilizing 
starand the time scale to collapse will be set by angular momentum transport 
processes that occur on a time scale that is substantially longer than the dynamical 
time scale 


Tns = 0.2 ms 



(4) 


where pnuc = 2.65 x 10^^ gcm“^ is the nuclear saturation density. Relativistic simula¬ 
tions of double neutron star merger^^^HUl indicate the central remnant is a HMN^ 
unless the initial ADM mass of the binary exceeds a threshold mass of 1.35 times 
the maximum mass of a cold, non-rotating neutron star. The discovery of two mas¬ 
sive neutron stars, J1614-2230 with M^s = 1-97 ± 0.04 MqI^ and PSR J0348-I-0342 
with Mns = 2.01 ± 0.04 Mq,I^ now place this limit to Mthresh > 2.7 Mq so that 
-depending on the distribution of neutron star masses that is realized in binary 
systems- a large fraction of the mergers may go through such a metastable phase. 
Taking the current mass estimates for binary neutron star^^ at face value (but 
keep in mind that some estimates have substantial errors) only three out of the 
10 DNS systems could avoid such a metastable phase and would instead directly 
collapse to a black hole. The delay time between merger and bh formation also has 
a crucial impact on both the triggering of short GRBs and -since it substantially 
impacts on neutrino-driven winds- on the merger nucleosynthesis. These issues are 
further discussed in Sec. andA long-lived HMNS could also produce additional 
radiative sigatures and contributions to high energy cosmic rays.l^^^t^ 

In most cases a CBM leads to the formation of a “thick” accretion disk where the 
scale height H is comparable to the radius RdisiQ The dynamical time scale of such 
a disk is roughly 


"^dyn,disk 


27r 

Wk 


0.01s 


M 

2.5Mq 


- 1/2 


77disk 

100 km / ’ 


while its viscous time scale is 


"Nisc ^ 0.3 S 



fldisk /2AMoy/' /TwyV 
100 kmy \ M J ^ 3 J ’ 


(5) 


( 6 ) 


where we have assumed that the viscosity can be parametrized as a Shakura- 
Sunyaev-type dissipation.!^ 


^In the following, we will use the abbreviation HMNS collectively for both “hypermassive neutron 
stars” that will finally collapse into a black hole and also for very massive, but stable, neutron 
stars. 

®We do not distinguish here between “disk” and “torus” and use both words synonymously. 
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3. Neutrino emission 

Like supernovae, compact binary mergers emit a large fraction of the released grav¬ 
itational binding energy 10^^ erg) in the form of neutrinos. This results pre¬ 
dominantly in neutrinos in the energy range of ~ 20 MeV. Their moderate energies 
together with the steep energy dependence of the interaction cross sections make 
them hard to detecl|^ But if the merger also accelerates material to relativistic 
speeds, as expected if they are the sources of short GRBs (Sec. [^, one may also 
expect neutrinos of substantially larger energy. The same shocks that are thought 
to accelerate the electrons responsible for the prompt gamma-ray emission should 
also produce relativistic protons, which can produce high-energy neutrinos.^^ilf^ 
Internal shocks that can produce the observed ^ MeV photons, should also be 
able to produce so-called prompt neutrino emission with energies in excess of 10^^ 
gy lUHlI] High-energy neutrinos may also be produced via neutron-rich outflowJUlIZS 
and reverse shocks.^ Such neutrino production mechanisms, however, rely on GRB 
outflows having a sizeable baryonic component. If instead the outflow should be 
essentially Poynting flux, large neutrino fluxes would be unexpected, unless they 
are produced efficiently in the forward shockEMZH 

Our main focus here, however, is the lower neutrino energy channel coming di¬ 
rectly from the merger remnant. Until merger, the neutron stars are still in cold 
/3-equilibrium, since tidal dissipation causes only a raise to moderate tempera¬ 
ture^^ (T ~ 10® K) and the tidal interaction has practically no time to change 
the neutron-to-proton ratio by weak interactions, see the inset in Fig. for an 
illustration. The /3-equilibrium condition on the chemical potentials is given a^ 
Jin = Jlp + fle + Jlv^ 1 where the bar indicates that the particle mass contributions 
are included. Since in old neutron stars no neutrinos are present (their diffusion 
time is of order seconds, see below) the neutrino chemical potential vanishes and 
the difference between the nucleon chemical potentials equals the electron chemi¬ 
cal potential which, in turn, determines the electron fraction T/. which has typical 
values below 0.1. After the merger, we can expect to have a fraction of the virial 
temperature, Tvir ~ 25 MeV (M/2.5 Mq) (100 km/R), i.e. temperatures substan¬ 
tially beyond the electron-positron pair production threshold (rUeC^ = 0.511 MeV). 
Therefore, positron captures n p + drive the electron fraction to higher 

values and yield the copious emission of electron-type anti-neutrinos which are the 
dominant neutrino species. This is different from the core-collapse SN case that 
yields neutrino luminosities in a similar regime, but dominated by the from the 
neutronization process p + e ^ n + In the hot, high-density interior of a HMNS 
one also expects heavy lepton neutrinos to be produced.^^SHUl Since they can only 


^Keep in mind that so far neutrinos have been detected only once for a supernova,ISSI SN1987A, 
and compact binary mergers occur ~ 10^ times less frequently. 

^We discuss here only the simplest case where neutrons, protons and electrons are present in the 
neutron star. When the Fermi energies become large enough to allow for new particle species to 
appear, the condition has to be modified accordingly, see for example GlendenningPH 
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interact via neutral current reactions, they escape easier from hotter and denser 
regions and therefore possess higher energies. Clearly, for the conditions prevailing 
in the remnant (p > gcm“^ in the HMNS and ^ 10^° gcm“^ in the disk) 
photons are completely trapped on the relevant time sacles and neutrinos are the 
only cooling agents. 

In the following, we will assume that a HMNS with temperature Thmns = 20 MeV, 
a mass M = 2.5 Mq and with a radius R = 20 km is present, so that its average 
density is p « 1.5 x 10^^ gcm“^(M/2.5 M©) {R/20 km)“^. We further assume 
for our estimates the presence of a surrounding accretion disk of Tdisk = 5 MeV, 
Mdisk = 0.2 M©, a typical radius i?disk = 100 km, a typical density pdisk = 10^^ 
gcm“^ and an aspect ratio H/R = 1/3, where H is the characteristic disk height. 
To obtain an order of magnitude estimate for the neutrino properties, we assume 
for simplicity that the main source of opacity is provided by neutrinos scattering 
off nucleon^ with a cross-section given bjESl cr = where the 

reference cross-section is ag = 1.76 x 10“"*^'^ cm^. Therefore, the mean free patlQin 
the HMNS is 


\HMNS 


1 

ntr 


1.8 m 


R 


20 km J \2.5 Mq 


M 




60 MeV 


and the corresponding neutrino optical depth is0 


^HMNS 


R 


;^HMNS 


lO"* 


R 


20 km 


M 


2.5 MqJ V60 MeV 


E„ 


(7) 

( 8 ) 


i.e. the HMNS is very opaque to its own neutrinos. The characteristic neutrino 
diffusion time scale is theiCHEZl 


^HMNS q.^HMNS ^ 
Miff ~ yr 


2.4 s 


R 

20 km 


-1 


M 

2.5 Mq 


E, y 

60 MeV ) 


(9) 


i.e. on the HMNS dynamical time scale, see Eq. the neutrinos are efficiently 
“trapped”. If we assume that a thermal energy of AEth ~ O.lGM^/i? can in prin¬ 
ciple be emitted, one expects a HMNS luminosity of 


^HMNS 


AEth 

_HMNS 

Miff 


3 X 10®2 erg/s 


M \ 
2.5 Mq ) 


E, 

60 MeV 


-2 


( 10 ) 


and, applying similar estimates as before, one finds the mean free path in the disk 


« 42km 


\ 10i^gcm“3 


I 15MeV 


( 11 ) 


^For Ue the opacity related to the absorption by neutrons is even larger, but still of the same order. 
^The expected neutrino energies are moderate multiples of the matter temperature, usually given by 
the ratio of two Fermi-integrals at the local chemical potential. For our scalings we use Ei, ^ Sk^T, 
i.e. 60 MeV for the HMNS and 15 MeV for the disk. 

JWe abbreviate both the optical depth and typical time scales with the symbol r. Since each time 
it has unique subscript this should no lead to confusion. 
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and therefore the optical depth is r ^ ~ 1. Thus, the typical escape time 

for a neutrino is substantially shorter than the dynamical and viscous time scales 
of the disk, see Eqs. ([^ and This short cooling time scale implies that neutrino 
emission can only be maintained through the constant supply of thermal energy 
via accretion. The gravitational energy gained by the accretion of the disk is then 
AE ~ GMMdisk/R and the expected neutrino luminosity becomes 

2^disk ^ 1 

2 Tvisc 

1053 ^ f M \ 5 / Mdisk \ / Q \ / Rdisk \ ^ f 100 km \ 5 /20jcmV 

s \2.5Mq) \0.2MqJ Vo.05^ \ J \ i?disk J \ R 1 

Clearly, these are only order of magnitude estimates, for more precise results nu¬ 
merical simulations are required. To date only a few implementations of neutrino 
physics in merger simulation codes exist It is encouraging that despite the 

steep temperature dependence of the neutrino emission rates (the electron capture 
energy emission rates, for example,are oc T®), the predicted neutrino luminosi- 
tie^HHHmiZl are rather robust, for a 2 x 1.4 Mq system around 2 x 10®® erg/s with 
average energies around 10, 15 and 20 MeV for Vf,, Ue and Vx, respectively, and, 
as expected from the simple considerations, with the luminosity being dominated 
by De- It seems in particular that implementation details play less of a role than 
the physical ingredients such as the treament of gravity or the equation of state. 
While the quoted results are based on relatively simple leakage schemes, more so¬ 
phisticated neutrino treatment^HEU yield actually similar results. It had however, 
been realized by Dessart et al.^ that nucleon-nucleon bremsstrahlung, which had 
been ignored in early implementations, is also an important source of heavy lepton 
neutrinos. 

More recently, neutrino emission has also been included in nsbh-merger simula¬ 
tions Ellllll Here the parameter range, especially the expected mass ratio range, is in 
principle even larger (and to date not yet well explored), but the first studied find, 
for Mbh = 7 Mq and large bh spin (obh > 0.7), also luminosities around ~ 2 x 10®® 
erg/s. Lower mass black holes are expected to yield hotter disks and therefore larger 
neutrino luminosities. Indeed, a simulation for a 5.6 Mq bh with Obh = 0.9 yield 
peak luminosities > 10®'^ erg/s. The bh mass distribution in nsbh binaries, however, 
is not well known and a recent studj®! based on X-ray binaries seems to suggest a 
peak around ~ 8 Mq . These results are consistent with studie^^ that find, from a 
binary evolution point of view, bh masses near ~ 10 Mq most likely. 

4. Compact binary mergers as cosmic factories of the heaviest 
elements 

The heaviest nuclei in the Universe form via neutron captures. In this process there 
is a competition between captures of new neutrons and /3—decays and the ratio of 
the two time scales can be used to distinguish between a “slow” (s-process) and 
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a “rapid neutron capture process” (r-process)piH^ Each process is responsible for 
about half of the elements heavier than iron. Although the physical mechanisms 
have been understood decades ago, there is still no concensus about the astrophysi- 
cal production site of the r-process. For many years supernovae have been considered 
as the most promising site, in particular the neutrino-driven wind that emerges after 
the formation of a hot, still deleptonizing proto-neutron star.l^SHinil Over the last 
decade, however, via increasingly sophisticated supernova simulation mo del^i^SHlfiSI 
the concensus has grown that the neutrino-driven winds do not provide the right 
physical conditions for the production of the third r-process peak (A ~ 195), but 
contributions to the less heavy region up to the second peak near A ~ 130 may 
be plausibly expected, see Arcones and ThielemamflSZl for an excellent review. A 
possible alternative supernovae site could be magnetohydrodynamically launched 
jets,112SI but these require a rather exreme combination of spin and magnetic field 
as initial conditions and at present it is not clear whether or at which rate such 
conditions are realized in nature. 

Neutron-rich matter ejected in a compact binary merger has been considered as an 
alternative formation channel of rapid neutron capture (”r-process”) elements since 
the 70iesl^51®li2£l Despite its physical plausibility, this channel was considered exotic 
for a long time and only the second-best model after core-collapse supernovae. This 
perception has only changed in the last one and a half decades through a number 
of studie^liiSKliZl that found r-process, up to the third r-process peak {A ^ 195), 
naturally occuring in compact binary merger ejecta. The ease with which r-process 
is produced by compact binary mergers together with the serious difficulties to pro¬ 
duce heavy r-process in core-collapse supernova^i^^Hinil lead to a shift in the 
general opinion towards compact binary mergers as the most likely production site 
of at least the heaviest, but possibly of even all of the r-process elements. An impor¬ 
tant issue that needs to be better understood, however, is whether/how r-process 
from compact binary mergers is consistent with galactochemical evolution. 


4.1. Mass loss chanels 

During the merger of a compact binary system mass is ejected into space via a 
number of channels. A fraction is released dynamically, i.e. via hydrodynamic inter¬ 
action and gravitational torques, and will subsequently be referred to as “dynamic 
ejecta”. Most CBMs also result in an accretion disk, and as such a disk evolves on 
a viscous time scale, a substantial fraction of its initial mass becomes unbound as 
a result of nuclear and viscous actionEZElHIll Another channel that has received 
a fair amount of attention recently are neutrino-driven winds jZHIZlir^l ^^ 

Highly magnetized neutron star-like merger remnants can also magnetically drive 
winds .E2IH1301 Also high-velocity ejecta coming from the interaction region between 
two neutron stars have been postulatedlHlJli23 All these channels differ in the 
amount of ejected mass, the distribution of electron fractions and entropies. There- 
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fore, they produce likely different element distributions and -if all produce electro¬ 
magnetic transients- these may be different for each channel. 

What is an ejecta amount that is interesting from a galactic r-process enrichment 
perspective? If TZi denotes the rate of event type i (averaged over the Galactic age) 
and ffii^c is the average ejected mass per event in a particular channel c of this event 
type, then the average Galactic r-process enrichment rate is 

•Adr.gai = TZ^rfij, (13) 

* c 

where Air gai is the r-process enrichment rate, averaged over the age of the Galaxy 
Tgai- In the following estimates we will conservatively restrict our considerations to 
nucleon numbers A > 140, since a large number of studies has shown that dynamic 
ejecta robustly produce such heavy r-process elements. It is, however, also a possi¬ 
bility that CBM produce all all r-process materialH^llill and it is also possible that 
event types different from CBMs contribute to the r-process production. This latter 
possibility will be ignored in our below estimates, therefore the quoted numbers 
should be interpreted as upper limits. 

If we assume that the solar abundance patterrP^ jg representative for the Milky 
Way, we can multiply the mass fraction of nucleP^ with A > 140, = 

2.6 X 10“®, with the baryonic mass of the Milky Way,li^ 6 x 10^° Mq, to find 
a Galactic r-process mass {A > 140) of « 1560 Mq and therefore the aver¬ 
age r-process enrichment rate {A > 140) is « 1.1 x 10“^ Mq year“^. Since 


Fe K 0.31 - 0.35 
s Kii 15 — 20kB/bar. 
V ^ 0.09 c 



0.10-0.18 c 


; central object: accretion disk 

dynamic ejecta HMNS or BH 


Fig. 3. A schemati c illust ration of the merger remnant. The numerical values are from recent 
numerical studies 
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nsbh mergers are estimated to be at least an order of magnitude rarer than nsns 
mergerswe only consider the latter for the following order of magnitude esti¬ 
mate Using’’plausible pessimistic” and ’’plausible optimistic” rate^^^as brackets, 
- 1 MWEG-i Myr-i and ~ 1000 MWEG"! Myr-\ where MWEG 

refers to “Milky Way-equivalent Galaxy”, we see that 

mnsns ~ M>l^°/TZ^sns ~ 10"^ ... 0.1 Mq (14) 

is an interesting ejecta amount from a chemical evolution perspective. If instead of 
only nucleon numbers A > 140 all r-process elements should be produced in GBM, 
the ejecta masses would need to be larger by a factor of « 13, if there are 

substantial contributions from other events the nsns yields need to be correspond¬ 
ingly lower. These estimates are what hydrodynamic simulation results need to be 
compared with. It further needs to be understood whether such masses/rates are 
consistent with the observed chemogalactic evolution. 


4.1.1. Dynamic ejecta 

In fact, essentially all recent binary merger calculations find dynamic ejecta masses 
in a range that is consistent with being a mjor source of cosmic r-process, see Table 
[l] For nsns systems, these dynamic ejecta can be further split into an “interaction 

Table 1. Comparison of the masses for the dynamic ejecta found in different numerical studies. While the 
numbers give a good impression about the expected mass range, it is worth keeping in mind that they are 
of limited comparability since the studies vary different parameters. CF: conformal flatness approximation, 
GR: General Relativity, PW: Paczynski-Wiita potential. 


authors 


reference dyn. ejecta [0.01 Mq] binary 


comment 


CF-approximation 

CF-approximation 


Bauswein et al. 2013 
Bauswein et al. 2014 
Deaton et al. 2013 
Foucart et al. 2013 
Foucart et al. 2014 
Hotokezaka et al. 2013 
Kyutoku et al. 2013 
Oechslin et al. 2007 
Rosswog et al. 2013 
Rosswog 2005 


1.7 X 10- 

2 ... 1.8 

nsns 

< 2 X 10- 

4 ... 9.6 

nsbh 

8 


nsbh 

1 .. 

5 

nsbh 

5 .. 

20 

nsbh 

0.01 .. 

. 1.4 

nsns 

1 .. 

7 

nsbh 

0.1 .. 

4.5 

nsns 

0.76 . 

.. 3.9 

nsns 

1 .. 

20 

nsbh 


Mbh = 5.6 Mq, Ubh = 0.9 
Mbh = 10 Mq, Ubh = 0.9 
Mbh — 7...10 Mq, Ubh — 0.7...0.9 
full GR 
full GR 

CF-approximation 

Newt. 

Mbh ^ 14 Mq, PW-potential 


component” that emerges from the interface between the two neutron stars and a 
“tidal component”f] that is launched by gravitational torques without experiencing 
noticeable shocks or shear flowsThe matter from the interaction region is 


’'Using simulation results on ejected mass,C22l the nsbh rate has recently been constrained to 
T^nsbh < 60 MWEG-i Myr-b 

’For an example see Fig. 2 in Korobkin et al. 2012.^^ 
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considerably hotter than the tidal component. Therefore, it may, via positron cap¬ 
tures, change its electron fraction Ye, while the tidal component is ejected at its 
low, initial /3-equilibrium value. The ratio of the mass in these components also 
depends on the dynamic evolution at merger and therefore on the treatment of 
gravity and the used equation of state. In Newtonian calculation^HEHEHl with the 
stiff Shen et al. EOSi^MUll the tidal component dominates while simulations using 
the Conformal Flatness approximation to GRpSl 

and simulations using dynamical, 
general-relativistic space-timeJ^Olare dominated by the interaction component. Due 
to their possibly different electron fractions this may have implications for nucle¬ 
osynthesis, see below. 

4.1.2. Disk dissolution 

When an accretion disk forms in the aftermath of a compact binary merger, it is 
characterized by p ~ 10^° gcm“^ , T ~ 10^° K and initial electron fractions close 
to the /3-equilibrium value, Yg < 0.1, of the disrupted neutron star, but increasing 
as the disk evolves viscously. Under such conditions the bulk of the disk consists of 
free neutrons and protons and e“'"e“-pairs. As the disk evolves viscously on a secular 
time scale, see Eq. ®, it expands and cools and at some point the free nucleons 
combine into light elements, predominantly alpha particle^^ which releases « 7 
MeV/nucleon. Together with vsicous dissipation, turbulent energy transport and 
energy deposition by neutrinos this can unbind a substantial fraction (up to 25 
%) of the original disk massIHZEIlHIllMIlIl How much matter becomes unbound by 
neutrino absorption is crucially dependent on the the presence of a central HMNS 
vs a black hole, in the latter case substantially less material is blown off from the 
H;.p[78l[79l[TT7P22l|T^ 


4.1.3. Neutrino-driven winds 

When two neutron stars merge, they release a substantial fraction of the gravi¬ 
tational binding energy in the form of neutrinos. As shown in Sec. the remnant 
neutrino luminosities are in the range of a few times 10®^ erg/s with typical neutrino 
energies of ~ 15 MeVEBEIESlllElll The gravitational binding energies per nucleon 
in the merged remnant, in comparison, are of order Egrav ~ ~35 MeV (Mco/2.5 
Mq) (100 km/r), where Mco is the mass of the central object. That means captur¬ 
ing a few neutrinos provides enough energy to potentially escape the gravitational 
attraction of the remnant. It had been realised early orEHEHElIl that the absorp¬ 
tion of neutrino energy could drive strong baryonic winds, similar to the case of 
a new-born proto-neutron starjUlUZl see Fig. for a simple sketch, but the wind 
properties (mass, electron fraction, entropy or geometry) are difficult to estimate 
without numerical simulations. 


•^Keep in mind that on the viscous time scale weak interactions may have substantially shifted 
the electron fraction Ye to larger values. 
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Early investigations used either order-of-magnitude approaches or parametrized 
models to explore the neutrino-winds from merger rernnantsE^mH E^ I ^^^ H^dHtl^ 
Due to the involved technical challenges, full-fledged numerical simulations of 
neutrino-driven winds have only been performed recently. The first neutrino- 
hydrodynamic study was performed in 2D by Dessart et al.,^ followed by re¬ 
cent studies by Fernandez and MetzgeiE^EdU and by Just et al.*JJ^ Just et al. 
explored the different ejecta channels using viscous hydrodynamics with Newtonian 
and pseudo-Newtonian gravity and they applied an energy-dependent two-moment 
closure scheme for the transport of electron neutrinos and anti-neutrinos. For black 
hole disk systems, they found neutrino interactions to help unbind disk material 
(~ 1% of the initial disk mass), however, to a much smaller extent than viscous 
action which is able to unbind up to 25% of the initial disk mass. This result is 
consistent with earlier studies in 2jjm 

that also found that neutrino heating in 
bh-disk systems unbinds only a small fraction of the disk material. 

To date, we are aware of only a single study in 3D, performed by Perego et al.lJ^ 
They explored the case where a HMNS is present and used the end points of 3D 
SPH simulations^ as initial condition. The further evolution was followed by means 
of an Eulerian hydrodynamics codd^^ augmented by a detailed, spectral leakage 
scheme that has been gauged at Boltzmann transport calculations of supernovae. 
They found a strong baryonic wind being blown out along the original binary ro¬ 
tation axis within ~ 100 ms, in qualitative agreement with the earlier 2D work 
by Dessart et al.l^ This wind unbinds at least 3.5 x 10“^ Mq and is therefore of 


relevance for the galactic chemical evolution, see Eq. (14). The weak interactions in 


the wind yield a broad distribution of electron fractions between 0.2 and 0.4, with 
tendentially more proton-rich outflow (Yg > 0.3) along the rotation axis while the 
equatorial outflow remained more neutron-rich (0.2 < Yg < 0.3). 


4.2. Nucleosynthesis in different channels 

All of the discussed mass loss channels are possible r-process sources. Most ad¬ 
vanced are to date the calculations for the dynamic ejecta, the other channels are 
still in early exploration stages. Already the early network calculations for the dy¬ 
namic ejecta of compact binary mergerd^J^J showed a very good agreement with the 
observed solar system abundance pattern, see Fig. left panel. Since no weak in¬ 
teractions were included in these calculations, the electron fraction Yg was varied in 
a low range considered appropriate for neutron stars. These results showed a very 
good agreement with the solar system pattern beyond A « 130 and in particular 
produced the heaviest elements up to the platinum peak without any fine tuning 
of the parameters, but hardly produced yields below A = 130. Such calculations 
have been more and more refined over recent yeardlMl^^^ll^^^lll^dlll^ and they all 
agree that the dynamic ejecta are a very promising site for (at least) the heaviest 
elements {A > 130). The abundance pattern has been showrPJ^to be virtually in¬ 
dependent of the astrophysical parameters of the merging binary system (masses. 
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mass ratio, whether the binary system consists of two neutron stars or a neutron 
star and a black hole) so that every CBM produces the same, very robust abun¬ 
dance pattern. This property is, of course, very interesting from a galactochemical 
point of view, since such robust abundance patterns are actually observed for the 
heaviest elements {Z > 56) in metal-poor stars.*^^ Moreover, studies that vary the 
neutron star EO^ come to the conclusion that the abundance is also rather sta¬ 
ble against such variations. A recent studjff^ has varied nuclear mass models and 
investigated the reason for the robustness. The authors suggest that a requirement 
for the robustness is that at the moment of freeze-out a much larger mass is in the 
fissioning region (nucleon number A > 250) than in the second r-process peak and 
above (120 < A < 180). 

While the results are -for fixed nuclear physics ingredients- very robust with respect 
to a change of the astrophysical parameters, they show a substantial sensitivity with 
respect to the variation of the nuclear physics input, they are particularly sensitive 
to the distribution of the fission product^HEHlimi which have a particular influence 
on the abundance pattern after the second peak. Although the abundance pattern 
globally agrees well with the solar system abundances, a close inspection shows in a 
number of studie^Hliii^lii^li^Ml^ZllIlll a slight, but systematic shift (with respect to 
the solar pattern) of the calculated results towards higher A. These are attributed 
to “late” neutrons, e.g. from /3-delayed neutron emission, but the agreement can 
be improvedl^li^ if ,d-decay rates are applied that are faster than those derived 
from the frequently used Finite Range Droplet ModeL^^ Such increased rates are 
actually consistent with recent experimental data for nuclei close to the r-process 
pattP^ and recent shell model calculations 

The nucleosynthesis in the other channels has so far been much less explored. In 



Massnumber A 


Fig. 4. Le ft panel: r-process abundance pattern found the dynamic ejecta from neutron star 
merger^i^ are shown as line, the solar system abundance pattern as crosses. Right panel: sketch 
of different event locations with respect to their host galaxy. 


both the neutrino-driven wind and the disk-dissolution channel matter stays much 
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longer, roughly a viscous time scale Eq. (|6]J^ at larger temperatures and in a neu¬ 
trino background field. Therefore, the weak interactions have time to substantially 
increase the electron fractions of the finally unbound material. The first hydrody¬ 
namic studies of these channel^^liill find that they produce less heavy r-process 
elements and with “non-robust” abundance patterns that depend on the exact de¬ 
tails on the merged system. Again, such a behaviour is consistent with the observa¬ 
tion of metal-poor starJi^ that show substantial deviations from the solar system 
pattern for .Z < 56 on a star-by-star basis. The wind- and disk-yields complement 
the dynamic ejecta at nucleon numbers below A « 130 and show broad consistency 
with the solar abundance pattern. 

Relativistic gravity also seems to have a substantial influence on the resulting nucle¬ 
osynthesis. This is mainly because neutron stars are more compact and the dynam¬ 
ics close to merger is impacted by GR effects. Overall, one expects a more violent 
dynamics, larger temperatures in the merged remnant and therefore an enhanced 
likelyhood to achieve a broader distribution of Yg values. A recent studj®^ includ¬ 
ing general relativity and a detailed neutrino treatment reports, for example, on a 
reproduction of the whole r-process range, however, without discussing in detail the 
channels and mechanisms. Closely related to the impact of GR is the equation of 
state and the GR-effects can be enhanced by particularly soft equations of state, 
while stiffer ones yield results that are closer to those employing Newtonian gravity. 


4.3. Galactic chemical evolution 


As outlined in Sec. 4.2 all recent studies based on hydrodynamic source plus nu¬ 
clear network simulations come to the conclusion that compact binary mergers are 
excellent candidates for a major r-process source, at the very least for the heav¬ 
iest nuclei (A > 130), but maybe even for all r-process nuclei. But there are, of 
course, also constraints from galactochemical evolution and the observed elemental 
scatter in metal-poor stars should be different depending on whether the majority 
of r-process comes from core-collapse supernovae or from compact binary mergers. 
For example, as discussed in Sec. |4.1| the typical ejecta mass of GBM of ^ 0.01 
Mq would -together with typical rate estimates- be able to account for most/all 
r-process material. If core-collapse supernovae, in contrast, would be main source 
and each supernovae would contribute, the average ejecta mass per event would 
only be ^ 10“^ Mq. Moreover, as will be discussed in the sGRB context, see Sec.[^ 
the large sensitivity of the gravitational wave inspiral times to the initial binary 
separation and eccentricity, see Eq. and Fig. naturally leads to a large spread 
in delay times after star formation. If endowed with substantial kick velocities, one 
expects a large variety of merger sites: some would merge close to their birth places 
in the galactic plane while others would merge in the outskirts of galaxies, see the 


^In cases where a stable neutron star results the time scales could be as long as the neutrino 
cooling time scales of several seconds, see Eq. ([9jl. 
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right panel of Fig. [^for a sketch. Such a distribution is indeed observed for short 
GRBs (also thought to result from CBM), they show typical projected offsets of 
~ 5 kpc, see Sec. |5.1| This also implies a large range of ambient matter densities 
encountered for different merger sites and therefore a different “stopping length” 
and volume with which the ejecta can mix. Thus, the galactochemical imprint of 
SN and CBM is likely quite different. 

Initial inhomogeneous chemical evolution studieJ^^ ruled out neutron star mergers 
as dominant source of r-process, mainly since -due to the low rates- the r-process 
enrichment would not be consistent with the observations at low metallicity and 
the scatter in [r-process/Fe] would be much too large. This is influenced on the one 
hand by the delay time due to CW inspiral, see Eq. ([^ and Fig.[^ and on the other 
hand by the question whether the neutron star ejecta mix with the chemical prod¬ 
ucts from the supernovae which formed the neutron stars in the first place. There 
are, however, indications that at least some of those supernovae that form neutron 
stars in surviving double neutron star systems are “non-standard”, they may form 
low-mass, low-kick neutron stars and eject only very small amounts of heavy ele¬ 
ments.More recent chemical evolution studie^^^ come to the conclusion that 
neutron star mergers could well be a major contributor to the galactic r-process in¬ 
ventory, although probably not the only one. Another chemical evolution studjES 
concludes that, except for the earliest evolutionary phases, CBM could well be the 
major production sites of r-process elements. 

Most recently, hydrodynamic studies of galactic enrichments have become available. 
In particular a high resolution cosmological hydrodynamics simulation called “Eris” 
has been used that is thought to closely represent the evolution of the Milky Way.l^^ 
The results of this simulation have been used in a post-processing step to follow 
the evolution of r-process enrichment by either type II supernovae or compact bi¬ 
nary mergers. Unlike in previous studies, the authors find that the nucleosynthetic 
products from compact binary mergers can be incorporated into stars of very low 
metallicity and at early times, even with a minimum delay time of 100 Myr and 
come to the conclusion that compact binary mergers could be the dominant source 
of r-process in the Calaxy. These results are supported by an independent hydro- 
dynamic studjff^ that uses cosmological zoom-in simulations of a Milky Way-mass 
galaxy from the Feedback In Realistic Environments project Also here, the au¬ 
thors conclude that the results are consistent with ns mergers being the source of 
most of the r-process nuclei in the Universe. 

Such an interpretation is also consistent with a recent studjff^of radioactive ^"^^Pu 
in Earth’s deep-see reservoirs. The authors find that the abundances are two or¬ 
ders of magnitude lower than what is expected from a continuous enrichment by 
supernovae. This points to a very rare actinide production site, for example a small 
subset of actinide-producing supernovae or a compact binary merger origin. 
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4.4. “Macronovae”: electromagnetic transients powered by 

radioactive decays from freshly syntesized heavy elements 

As outlined above, compact binary mergers eject, via several channels, initially ex¬ 
tremely neutron-rich matter. The first and best-studied channel are the dynamic 
ejecta. Li and PaczynskP^ were the first to realize that if indeed rapid neutron 
capture occurs in the ejected materia£3i22l[III] then the radioactivity in the ejecta 
should cause an electromagnetic transient following a merger. Such transients are 
often referred to as “kilonovae’’^^ or “macronovae” But since only little mass 
it ejected (ruej ^ 0.01 M 0 ) at large velocities (uej ^ 0.1 c), the expected light curve 
evolution is faster and the peak luminosity is lower than in the type la case. The bulk 
of the energy release due to r-process occurs on a time scale of ~ secondliiliii^li^liZSl 
when the matter density is still so large that radiation is completely trapped and 
the radioactive heating is used up for accelerating the adiabatic expansion of the 
ejectaj^ It is only when the expansion time scale becomes comparable to the radia¬ 
tive diffusion time scale that substantial electromagnetic emission occurs.^i^ The 
diffusion time scale and therefore the time to transparency depend on how opaque 
the material is with respect to its radiation content, but unfortunately the opacities 
are to date maybe the largest uncertainty in modeling macronovae. 

For simplicity, the first study of Li and Paczynski (and all subsequent studies until 
recently) assumed for a lack of better knowledge that the relevant opacities would 
be of order k ~ 0.1 cm^/g, similar to the line opacities of iron-group elements.li^ 
The composition of the dynamic ejecta, however, is very different from any other 
known cosmic explosion and in particular very different from any type of supernova. 
While the latter produce elements up to the iron group near Z = 26, the dynamic 
ejecta of neutron star mergers consist, due to the extreme initial neutron-richness 
of the decompressed r-process material, almost entirely of r-process elements up to 
the third peak near Z « 90, e.g. Fig. To date, not much is known about the 
relevant opacities of such material, but it was recently realized by Kasen and col- 
laborator^i^liHl that the opacities have previously been underestimated by orders 
of magnitude. The literally millions of lines, Doppler-broadened by the differential 
velocities of the ejecta, contribute a pseudo-continuum of bound-bound opacity. 
Photons trying to escape the ejecta come into resonance with multiple transitions 
one by one, resulting in a very large opacity. The bulk of the opacity is thereby 
provided by ions with a particularly complex structure. In particular lanthanides 
(58 < Z < 72) have -due to the complicated structure of their valence f shells- been 
found to be major sources of opacity. While detailed knowledge of opacities is still 
lacking, the currently existing calculations indicate opacities that are ^ two orders 
of magnitude larger than what was previously assumed ,k ^ 10 cm^/g, even if 
only small amounts of lanthanides are present. As a result, radiation is much longer 
trapped and matter becomes transparent only much later, at lower temperatures 


For the nuclear heating history see, for example, Korobkin et al.,E^ their Fig.7. 
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and luminosities and peaking in the near-IR band rather than, as initially found, in 
the optical. Using such increased opacities, recent studieJ^^^HIlSI found the resulting 
“macronovae” should peak after typically a week rather than the original estimate 
of half a day. 

Simple order of magnitude estimate£13Ell[Illl can be obtained in the following way. 
The diffusion time scale for the ejecta with radius R is approximatelj]^ 


^difT 


rrie^K 

AttcR 


(15) 


The peak emission is expectecP^ when this is equal to the expansion time scale 
Texp = so at a radius 
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If we assume that the late-time radioactive energy release rate can be approximated 
as a power lawpZlEMiZll e - eo(iAo)-“ (a = 1.3), we can estimate the peak 
bolometric luminosity as 


-'peak 
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where we have used the numerical values as determined in Korobkin et al. (2012). 
Using the Stefan-Boltzmann law, L = Airas bT‘^■, one can obtain a rough estimate 
for the effective temperature at peak as 


Tpfak = 2200K 


^10 cm^/g^ 
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( 20 ) 


These simple estimates also illustrate the dramatic effect that the changed opacities 
have: instead of a transient that peaks in the blue part of the optical spectrum after 
~ 1 day (for iron-type opacities), one now expects a peak after ~ 1 week in the 
near-lR. 

Early optical searches for macronova emission were carried out for several short 
GfoBsEsmUIl but none of them provided convincing evidence for a macronova. The 
first serious evidence came from a search in the near-IR band in the aftermath of 
GRB 130603B using a combination of ground-based observations at < 2 days and 
HST observations at about 9 and 30 days after the burst EHUEUI The difference be¬ 
tween the observations revealed the presence of a near-IR point source at 9 days 


PSee, e.g., Piran et al. (2013).C^ 
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after the GRB, this source had disappeared by the time of the second observation. 
These observation lead both group^HlElll to conclude that the most natural ex¬ 
planation was a radioactively powered transient due to decaying, freshly produced 
r-process elements, in other words by a macronova. 

While overall compatible with the expectations for a macronova caused by the 
extremely neutron-rich ejecta, the inferred value^MEHlElll for the ejecta masses 
(ruej « 0.03 — 0.08 M©) and velocitieJ^^ (uej « 0.2c) are in a plausible, though 
slightly uncomfortably large regime for a near-equal mass nsns merger which is con¬ 
sidered the most likely event. A natural way out would be a CBM with a mass ratio 
substantially different from unity, either nsns or nsbh, which would explain both 
larger ejecta masses and velodtieslii^llli^ Such systems, however, are expected to 
occur less frequently than g « 1 nsns-mergers. While r-process powered macronova 
are a natural interpretation, there is still room for alternative model£13I13 and 
further macronova detections are eagerly awaited to settle the case. 


5. Compact binary mergers as central engines of short gamma-ray 
bursts 

Gamma-ray bursts are intense flashes of soft gamma-rays that are detectable ap¬ 
proximately once per day and that reach Earth from random directions. There were 
early hint a structured duration distribution, but this became only firmly 

established through the work of Kouveliotou et who found a minimum in the 

duration distribution around ~ 2 s. They also realized that bursts below this time 
scale are consistently harder in their spectra than the bursts above it. This lead to 
the establishment of two classes of bursts, ’’short, hard bursts” (sGRBs; r ^ 0.3 
s) and ’’long, soft bursts” (IGRBs; r ^ 30 s)j^ Their spectra in the conventional 
ly — vFjj coordinates are usually fitted by a broken power law with a pronounced 
peak around ^ 400 KeV energies (for sGRBs; ^ 200 keV for IGRBs) 

Compact binary mergers consisting of two neutron stars or a neutron star and a 
black hole have been suggested as one of the first progenitor models ji2£ll22IiEII] ac¬ 
tually long before the bimodality of the GRB distribution,l22iI their cosmological 
origiilHIllHiHlIZi^ (and therefore the distance and energy scales) and their types of 
host galaxie^^UHSm -^^ere firmly established. The CBM model has survived the con¬ 
frontation with three decades of new observations and many of the arguments that 
were brought forward in the very early papers are still considered valid today. 

The discussion here will be centered around the role of compact binary mergers, for 
more exhaustive, general discussions of GRBs and their properties, we refer to the 
excellent review^^UHSH] exist on the topic. 


‘iThis classification has, however, been subject to criticismll2ll22ll 





















February 12, 2015 1:35 


ms 


22 S. Rosswog 

5.1. Confronting the CBM model with observations 

We will here summarise to which extent CBMs explain the observed properties of 
sGRBs and where tensions between the model and observations exist. Alternative 
suggestions are briefly discussed in Sec. |5.5[ 

Energy requirements 

Although often paraphrased as ’’the biggest explosions in the Universe” the en¬ 
ergy requirements of sGRBs are large, but still moderate compared to the rest 
mass energy of a solar mass, 1 MqC^ « 1.8 x 10^"^ erg. The observed, ’’isotropic” 
7 -energies (i.e. assuming that the radiation is emitted isotropically) are in the 
range of = 10 ^®... 10 ®^ erg, 1225 H 22 l] -vvhich correspond to ’’true” energies, if 

the radiation is emitted into a solid angle with half-opening angle 9, of E^ = 
10^®...10®° erg (0/8°)^. In fact, for the small subset of sGRBs where opening angle 
information is available the corrected values for both 7 —ray and kinetic energies 
indicate E^ ~ Akin ~ erg.l22SI The main energy reservoir that can be tapped in 
a GBM is the released gravitational energy which is Agrav ~ where Mtot 

its the total binary mass and R a separation not much larger than the Schwarzschild 
radius Rg = 2GMtot/c^ = 9 km (Mtot/SM©). So that, in principle, an energy ap¬ 
proaching ~ MtotC^/2, corresponding to several times 10®® erg, is available. There¬ 
fore, the energy requirements are not a serious challenge, even if only moderate 
efficiencies should be applicable for transforming released energy into 7 -rays. 

Variability time scale 

The energy output of sGRBs can vary substantially on a time scale of a few 
msEHHHZI This is the time scale that is naturally expected from a compact binary 
system. Both from the dynamical time scale of a neutron star, see Eq. Q, and from 
the orbital period at the innermost, stable circular orbit (ISCO) of a stellar-mass 
black hole, see Eq. ([^, one would naturally expect variations of order milliseconds. 

Duration 

It is commonly expected that an accretion disk is needed to transform the released 
gravitational binding energy into ultra-relativistic outflow and finally into the ob¬ 
served radiation. If one assumes that such extremely rapidly accreting systems can 
be described by the simple thick-disk estimate Eq. (§, there is good agreement with 
the typical duration of sGRBs,f 2 lII ^ o.3 s. 

Host galaxies/environment 

The host galaxies of short bursts are systematically different from those of long 
bursts. While the latter occur in unusually bright star-forming regionsthe 
former occur in both early- 20%) and late-type 80%) galaxies ESIEHII This 
points to a generally older stellar population and a substantial spread of ages. This 
is actually expected for compact binary mergers and has in fact been a prediction 
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of the merger modelP^ As discussed in Sec. the inspiral time due to the emission 
of gravitational waves is a sensitive function of the initial orbital period P and the 
eccentricity e, see Eq. Q and Fig. Since the initial values of both quantities are 
set by evolutionary processes of the stellar binary progenitor, one expects indeed 
a broad spread of inspiral times. Due to kicks that the neutron stars receive at 
birth large initial eccentricities may be commonj^If the binary as a whole receives 
systemic kick, substantial distances can be travelled during the inspiral time and 
one expects a broad distribution of merger offsets with respect to their host galax¬ 
ies. The distributions obtained from binary evolution calculation^HlElIElIllMllagree 
remarkably well with those observed for sGRBsEmEMElll The observed projected 
offsets range from 0.5 to 75 kpc with a median value of « 5 kpc. These results 
are interpreted as projected kick velocitie^^H! from Vkick ~ 20 to 140 km s“^ (a 
median value of ~ 60 km s“^) and consistent with the binary population synthesis 
results (ukick ~ 5 to 500 km s“^) for nsns mergers |2351t238| Also, the sGRB red- 
shifts expected from CBM models are in agreement with those 
0.1<z<1.5). 


Beaming and event rates 

The angular structure of the outflow can be inferred from so-called achromatic 
“jet-breaks” \23MM11 the light curve decay begins to steepen roughly simultaneously 
in all wavelength bands. This is caused by a combination of “relativistic beam- 
and, in addition, by the sideway expansion of the jetted out- 


If the outflow of sGRBs is collimatecQ into a half-opening angle 0, the emission 
from each emitting patch is “beamed” into a forward cond^iS) .^^ith half-opening an¬ 
gle 1 / 7 , where 7 is the local Lorentz factor. As a consequence, only observers that 
happen to be in the beam cone can observe the emission and for them it looks like a 
spherical flow with the properties of the local patch. While the jet spreads sideways 
and the outflow becomes decelerated by the interaction with the ambient medium, 
the relativistic beaming angle increases, and once similar to the true geometric col- 
limation angle, I /7 ~ 0 , the jetted nature of the outflow becomes noticeable and a 
jet break in the lightcurve indicates that so far only a small patch had been visible. 
If redshift, kinetic energy and ambient matter density are known, one can infer the 
geometric jet opening 9 from the time of the jet breal|^ Obviously, the value of this 


"^Note, however, that some systems may have formed via different mechanisms. In the case of 
the double pulsar system, PSR J0737-3039A/B, the low space velocity, the comparatively low 
mass of pulsar B, the low orbital eccentricity (e = 0.09) and the location only ~ 50 pc from the 
Galactic midplane are all consistent with pulsar B hav ing been formed by a mechanism different 
from the usually assumed core-collapse of a helium starll£IL^£J In fact, progenitor masses possibly 
substantially below 2 M 0 are kinematically favored. 

®We refer to “collimated” to describe the geometry of the outflow and to “beamed” for the special- 

relativistic aberration effect. _ 

^See, for example, Eq.(33) in Nakar (2007). 
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angle has severe implications for both the overall energy budget, see above, and the 
true event rate: for small opening angles, most bursts are unobservable and the true 
rate is larger by a “beaming factor” fh = Att/IWI « 2 / 0 ^, where AS! is the solid 
angle of the beam. 

So far, it has been difficult to observe jet breaks for sGRBs, mainly because of their 
weak afterglows. There is only a handful of cases for jet breaks and they have values 
of ~ 10°|2i51l2illin rough agreement with the theoretically expected values 
This would imply that only roughly one out of 70 sGRBs is detectable and it would 
translation! into a true sGRB rate of « 20 yr“^ within 200 Mpc, roughly the dis¬ 
tance accessible to Advanced LIGO/VIRGO. 

’’Macronova” emission 

“Macronovae”, see Sec. |4.4[ have actually been a prediction of the compact binary 
merger model. The detection of the first the macronova event in the aftermath of 
GRB1306033i2^li^ is broadly consistent with the expectations for radioactively 
powered transients that result from the extremely neutron-rich dynamic ejecta that 
produce very heavy nuclei {A > 130). If the current ideas about other mass loss 
channels, see Sec. |4.1[ are correct, there should be additional transients with differ¬ 
ent properties, also powered by radioactivity. 

Lorentz factors 

GRBs, both long and short ones, involve highly relativistic bulk flows. This was 
suggested as a solution of the so-called “compactness problem”l 21 £] 12 ^ and later also 
confirmed by direct observations.!^^ The clue to the compactness problem is to 
understand, how a luminous source can vary on time scales of order milliseconds 
and still emit optically thin radiation. Assume a (slowly moving) sphere of radius 
R that emits radiation and changes substantially on a time scale St. If it would 
“switch off” immediately, it would take at least the light travel time difference 
between closest and the farthest visible point of the sphere, R/c, to conveigh this 
information to an observer. Turned around, (assuming a non-relativistic source) an 
observer would conclude that the source must be of a size R < cSt cz 300 km(5t/I 
ms). Now inserting typical numbers of GRBs, such a size would imply enormous 
optical depths (r ^ 10 ^^) for photons with respect to the at photon energies of ^ 
MeV copiously present electron-positron pairs ,*2^ in stark contrast to the observed, 
optically thin radiatiorj^ 

If the source is instead moving relativistically to the observer, the source can be 
larger by a factor of ^ 7 ^ (because the observer can, due to relativistic beaming, 
only see a small patch of the source, but the real source size is much larger) and 
also the photon energy in the emission frame is lower than the observed one. The 
requirement r < 1 can then be used to place lower limits on the Lorenz factor, 
estimates typically favour 10^ < 7 < 10^, see e.g. Lithwick & SarP^for a detailed 


This argument had originally been used to place an upper distance limit on GRB sources. 
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analysis. 

A second line of argument comes from the onset of the afterglow. Once the blast 
wave has deposited most of its energy in an ambient medium, it assumes a self¬ 
similar profil^^m and X-ray and optical emissions decay as powerlaws. During this 
powerlaw phase the evolution of the Lorentz factor (only weakly depending on the 
kinetic energy of the blast and the ambient medium density) is known, and the time 
of the onset of this powerlaw behaviour provides a lower limit on the Lorentz factor. 
The presence of such large Lorentz-factors, is actually a very non-trivial constraint 
on GRB models. To accelerate to an assymptotic Lorentz factor of 7 , a blast of 
energy if, cannot contain more than 

mtar = 2L = 5.5 X 10-®Mo (, (21) 

1C V7100/ 

where we refer to a quantity A„ as A/10” in cgs-units. It is at least not obvious, 
how a small fraction of the mass can receive such a disproportionate share of the 
released energy. How are the deposition of mass and energy so cleanly separated? 

”Late-time activity” 

There are also a number events that show “late-time activity” on a time scale much 
longer than what is naively expected from the compact binary merger model, see 
Sec. To date this activity is still only incompletely understood. 

A subset of ^ 20 % of sGRB^HlE^ shows after the first, “prompt” spike somewhat 
softer, extended 7 -ray emission that lasts for ^ 10 - 100 s and is sometimes delayed 
in its onset. This emission component was first seen in a stacking analysis of short 
BATSE burst^^m and subsequently observed in a number of other burstsliH2ll2^^1t2^ 
In some cases, the fluence in this “extended emission” can exceed the one in the 
prompt spike, in the extreme example of GRB080503 the extended emission contin¬ 
ues up to ^ 200 s and dominates the fluence of the initial spike by a factor of 32 . 1^201 
Another type of late-time activity are X-ray flares following the prompt 7 -ray emis¬ 
sion which have actually been observed for both long an short bursts.l^snHUH 
some cases, such flares have been seen many hours after the burst, e.g. GRB050724 
showed a significant X-ray flare at ~ 14 hours after the burst, for GRB 130603B 
excess X-ray emission has been detected more than a day after the main burst. For 
sGRBs such flares occur for different types of GRB host galaxies and both for cases 
with extended emission and without.^^^ 

Such time scales that exceed both the dynamical and the viscous time scales by 
many orders of magnitude are clearly uncomfortable for the compact binary merger 
model and had not been expected prior to detection. Given that main power source 
of GRBs is accretion it is natural to start from the hypothesis that it may also power 
this late-time activity. The dissipation parameter a in Eq. is not well known, but 
it seems unlikely to have a value that maintains accretion for much longer than 
seconds. Gompact binary mergers (in particular those with mass ratios deviating 
from unity) also eject a few percent of a solar mass into nearly unbound orbits. 
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and this material while providing an energy reservoir > erg provides a time 
scale, Tfb, that can easily exceed both the dynamical and the viscous time scale 
by orders of magnitude I21311Z1I1111 However, at least in its simplest form, it may 
be challenged to release large amounts of energy at very late times (s> seconds). 
Dynamical collisions between compact objects, see Sec. |5.4| carry some promise in 
that respecll53lH2iJE®£lEZ2! since each of the several close encounters produces a tidal 
tail that serves as a mass reservoir far from the engine. Most likely, however, they 
are too rare to explain the common late-time activity phenomena. While the initial, 
purely ballistic fallback ideas may have been too simple to explain the observations, 
it has become clearer over the last decade that compact binary mergers eject mass 
via several channels, see Sec.|^ some driven by neutrino- and/or nuclear processes, 
and the interaction of these different mass loss channels may yield a much more com¬ 
plicated picture than just the ballistic fallback. In that sense, there is some promise 
that also the late-time activity can be reconciled with compact binary mergers, but 
this question is certainly far from being understood at the moment. 


5 . 2 . How to launch ultra-relativistic outflow 


Models for the ultra-relativistic outflow must -apart from questions of collimation 
and stability- first and foremost explain how some fraction of matter can acquire 
such a disproportionate share of the energy. The very large energy to rest mass ra¬ 


tio required to produce the inferred ultra-relativistic outflow, see Eq. (21), is a very 


non-trivial requirement and to date it is still a matter of debate how this is achieved 
in GRB. One of the mechanisms suggestedI^112£l[I?^^^^IEIUt2I£l to achieve this is the 
annihilation of neutrino anti-neutrino pairs, RiVi —)■ e“e“'". The annihilation cross- 
section is very small, but the neutrino luminosities are huge, see Sec. and if a 
small fraction of the neutrino energy can be deposited in a baryon free region the 
typical energy of a short GRB, ^ 10^® erg, see Sec. ^ could plausibly be reached. 
Due to the dependence of the annihilation rate on a factor ^ = 1 — cos 0, where 
0 is the collision angle between neutrino and anti-neutrinc|^ headon collisions are 
favoured. Therefore, spacetime curvature or thick disk geometries can enhance the 
annihilation rate. A necessary requirement for neutrinos from an accretion disk is 
that they can escape on short enough time scale so that they are not just advected 
into the bh. But as estimated in Sec. and further detailed by numerical simula¬ 
tions, the neutrino escape time scales are substantially shorter than the dynamic 
disk time scales. 

Large neutrino luminosities, however, can also lead to strong neutrino-driven winds 
that can make it impossible to reach the required energy to rest mass ratio. This 
effect is particularly pronounced for remnants that contain a central HMNSEHEU 
The situation is less critical for black hole disk systems, for which numerical simu- 


''In the explicit expression|^2IlZZl there are actually two such terms, one is proportional to p, the 
other to g?. 
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latioii^IZElll find much weaker neutrino-driven winds. A recent studjEzH came to 
the conclusion that the HMNS must collapse within ^ 100 ms, otherwise baryonic 
pollution precludes the emergence of ultra-relativistic outflow. 

Magnetic mechanisms are another broad class of plausible mechanisms. They are, for 
example, required to extract the energy of a spinning black hole via the Blandford- 
Znajek mechanism,or they can more directly be transformed in outflow by buoy¬ 
ancy instabilities with subsequent reconnection eventor as some form of 
extreme pulsaiEiaMHlll where the rotational energy is tapped. All these processes 
require very strong, near-equipartition magnetic field strengths to be effective. This 
is, however, likely to be achieved, at least for the short relevant time scale, since 
neutron stars are from the beginning endowed with strong magnetic fields and the 
merger process offers ample possiblities to amplify themllMi22IIH®lHlIZl 
Of course, it is certainly a possibility that the diversity of observed bursts is related 
to several of these mechanisms work in concert. 

5.3. Double neutron star vs neutron star black hole binaries 



Fig. 5. Left panel: illustration of the mass transfer separation a^T with respect to horizon and 
innermost stable circular orbits based on Newtonian estimates (for fiducial neutron st ar p roperties, 
Mns = 1-4 Mq and i2ns = 12 km). Right panel: results of an analytical modeP^ fitted to 
numerical simulations for the mass that remains at late times outside the hole. The remaining 
mass depends apart the black hole mass M^h and spin Xbh (= <^bh ia our notation) also on the 
neutron star c omp actness Cns- Results for more/less compact neutron stars can be found in the 
original paper.^^^ 


Traditionally, nsns and nsbh systems have been considered as “standard GRB” 
progenitor, often without much distinction between them, since they were thought 
to both lead to the most likely engine, a stellar mass bh surrounded by an accretion 
disk. In recent years, however, this picture has become more differentiated. On the 
one hand it has become clear that bh formation in the nsns case does not have 
to happen early on and may actually in some cases not happen at all, see the 
discussion in Sec. The presence of a HMNS, in turn, implies strong neutrino- 
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driven winds (Sec. 4.1.31, which, in turn, pose a potential threat to the ability of 
launching ultra-relativistic outflow. The nsbh case supposedly produces a baryon- 
cleaner environment (mainly due to the absence of a HMNS) and it is more likely to 
show imprints from jet precession.l^^ But not every combination of bh masses and 
spins is actually able to form a substantial accretion torus around the hole. This is 
mainly because the separation where the neutron star starts transferring mass/is 
disrupted grows with a lower power of the bh mass Mbh than the innermost stable 
circular orbit (ISCO). Therefore for large enough bh masses the neutron star is 
swallowed before being disrupted so that no torus can form. 

The final answer requires 3D numerical simulations with the relevant physics, but a 
qualitative idea can still be gained from simple estimates. Mass transfer is expected 
to set in when the Roche volume becomes comparable to the volume of the neutron 
star. By applying Paczynski’s estimate for the Roche lobe radiu^^HI and equating 
it with the ns radius, one finds that the onset of mass transfer (which we use here 
as a proxy for the tidal disruption radius) can be expected near a separation of 

V3 / p X 


OMT = 2.17i?n 


1 + 9 


26 km 


12 km 


1 + 9 


( 22 ) 


Since omt grows, in the limit where the mass ratio q = Mns/Mbh 1, only pro- 
portional to , but the ISCO and the event horizon grow oc Mbh, the onset of 
mass transfer/disruption can take place inside the ISCO for large bh masses. At 
the very high end of bh masses, the neutron star is swallowed as whole without be¬ 
ing disrupted at all. A qualitative illustration (for fiducial neutron star properties, 
Mns = 1-4 Mq and i?ns = 12 km) is shown in Fig.[^ left panel. Roughly, already for 
black holes near Afbh ~ 8 Mq the mass transfer/disruption occurs near the ISCO 
which makes it potentially difficult to form a massive torus from ns debris. So, low 
mass black holes are clearly preferred as CRB engines. Numerical simulation^^lll 
have shown, however, that even if the disruption occurs deep inside the ISCO this 
does not necessarily mean that all the matter is doomed to fall straight into the 
hole and a torus can still possibly form. 

The bh spin also plays a crucial role for the question of disk formation, since it 
determines the location of the ISCO, for maximally spinning Kerr bhs, for example, 
it is located at only Risco = GMbh/c^. This, of course, has a serious impact on the 
dynamics and in particular on the resulting tori. This is illustrated in Fig. where 
two nsbh binary systems are simulated (using the SPH method) that are identical 
(initially circular radii with Tc = 8GMui/c^ of a 1.3 Mq ns around a spatially fixed 
Kerr bh) up to the bh spin parameter. In the first case of a non-rotating bh (upper 
row, abh = 0) the ns is swallowed completely within ~ 13 ms, while for the rapidly 
spinning case {ubh = 0.9) a massive torus forms. 

Fully, relativistic simulations are computationally very expensive, but it is possible 
to devise simple analytical model^^SS] jj^^ss remaining outside the blj^ that 


'^This refers not only to disk mass, but also contains in addition dynamic ejecta and tidal tail. 
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have a physically intuitive functional form and whose parameters can be fitted to 
the results of relativistic nsbh simulations. Once fitted, the results depend on the 
black hole mass Mbh, its dimensionless spin parameter obh and the neutron star 
compactness Cns- The result for the probably most likely neutron star compactness 
is shown in Fig.[^ right panel (their quantity xbh is what we refer to as Obh)- The 
results for more extreme compactness values can be found in the original paper.^^lll 
According to this model, a 10 Mq bh needs to have at least a spin of Ubh ~ 0.9 to 
form a disk in the disruption of a “standard” neutron star, for larger/smaller stars 
(13.5 km/9.5 km) spin values of 0.84/0.98 are needed. 

When discussing disk formation in a GRB context, it is worth keeping in mind that 
even seemingly small disk masses allow, at least in principle, for the extraction of 
energies, 

Aextr ~ 1.8 X 10^1 erg (—) , (23) 




Fig. 6. Impact of black hole spin on disk formation. The panels shown in the upper row refer to 
the merger of a neutron star with a Schwarzschild bh (flbh = 0)? while the lower row shows the 
merger with a rapidly spinning Kerr bh (ath = 0.9). In both cases a 1.3 Mq neutron star is set on 
a circular orbit with Vc = 8GMbh/c^- Due to the ISCO lying at a larger radius, for the Ubh = 0 
case the black is swallowed almost completely (only very few SPH particles are left at t = 13 ms), 
while for the highly spinning case a massive accretion disk forms (last panel, second row). Note 
that the panels in the two rows show different scales. 
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that are large enough to accommodate the isotropic gamma-ray energies, E^^iso ~ 
IqSO g]-g^ have been inferred for short burstsP^^ Collimation into half-opening 
angle 0 further reduces the requirement by 0^/2. 

If the bh mass distribution in nsbh binary systems is indeed peakedlHlIini around 
8 — 10 Mq then the large spins required to form a disk (obh > 0.8) may seriously 
constrain the parameter space. Recent spin measurements in X-ray binarie^^MI 
in four out of ten cases spins > 0.8, thus giving some hope that such high spins may 
possibly be realized also in nsbh binaries. A recent studjPHsI finds that for most of 
their Monte Carlo models between 10 and 30% of the nsbh mergers could be able 
to launch a sGRB, but the fraction could essentially vanish for bottom-heavy bh 
spin distributions. One should therefore be prepared that the contribution of nsbh 
mergers to the observed sGRB rate could actually be small. 


5.4. Gravitational wave-driven binary mergers vs. dynamic 
collisions 

Traditionally, mostly gravitational wave-driven binary systems such as the Hulse- 
Taylor pulsaiEElHl have been discussed as GRB engines, dynamical collisions among 
neutron stars and black holes in stellar systems with large number densities such 
as Globular Clusters have been considered a very unlikely event More recently, 
however, such collisions have received a fair amount of attentionEIll2^llZ3E£M222l 
Unfortunately, their rates are at least as difficult to estimate as those of GW-driven 
mergers. 

Collisions differ from gravitational wave-driven mergers in a number of ways. For 
example, since gravitational wave emission of eccentric binaries efficiently removes 
angular momentum in comparison to energy, primordial binaries will have radiated 
away their eccentricity and will finally merge from a nearly circular orbit. On the 
contrary, binaries that have formed dynamically, say in a globular cluster, start from 
a small orbital separation, but with large eccentricities and may not have had the 
time to circularize until merger. This leads to pronouncedly different gravitational 
wave signatures, “chirping” signals of increasing frequency and amplitude for merg¬ 
ers and initially well-separated, repeated GW bursts that continue from minutes 
to days in the case of collisions. Moreover, compact binaries are strongly gravita¬ 
tionally bound at the onset of the dynamical merger phase while collisions have 
total orbital energies close to zero and need to shed energy and angular momentum 
via GW emission and/or through mass shedding episodes in order to form a single 
remnant. Due to the strong dependence on the impact parameter and the lack of 
strong constraints on it, one therefore expects a much larger variety of dynamical 
behavior for collisions than for mergers. 

A study that discussed various aspects of such collisions in detaiP^SI came to the 
conclusion that such encounters could possibly produce an interesting contribu¬ 
tion to the observed GRB rate. Indeed, such collisions typically yield large accre¬ 
tion tori masses, peak neutrino luminosities^ in excess of 10®^ erg/s, but in some 
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cases substantially more, and ample possibilities to amplify initial seed magnetic 
fieldsPerhaps most interestingly, since they typically suffer several 
close encounters before they can merge into a single object, and since each close 
encounter results in tidal tail, they produce a more complicated environment of the 
central merger remnant. This stores in particular subtantial amounts of mass far 
from the remnant which only falls back on longer time scales and is maybe related 
to the late-time activity of some sGRBs, see Sec. |5.1[ 

Collisions, however, also share a possible caveat with the merger cases: at least in 
the nsns case they likely produce very strong neutrino-driven winds that could pre¬ 
vent triggering a short GRB. More recently, it was concluded that collisions cannot 
be very frequent, since they eject very large amounts of neutron-rich matter. The 
dominant nsbh collision channej^ ejects typically 0.15 Mq per event which 
-from cosmic nucleosynthesis arguments- restricts the (nsns -I- nsbh) collision rate 
to a maximum of « 10% of the estimated nsns merger rate.l^The rare nsbh collision 
cases, however, would result in particularly bright “macronova” emission with 


Lpeak«10«erg/s (^) 


f /10 cm^/g 


nie 


0.15Mq 

peaking after approximately 13 days, see Eqs. ([T7| and (19). 


(24) 


5.5. Alternative models 

As an alternative to the common bh-plus-disk systems, magnetars have 
been suggested as central engines of sGRBs, both for long and for short 
bursts 111]119 iJ128ij[282]|302j| 3^ Magnetar-like objects could be formed under a variety 
of circumstances. Giant flares from “young”, extragalactic magnetars with long re¬ 
currence times could plausibly be responsible for a fraction of the events observed 
as sGRBs.l^2lll308| A progenitor population of young magnetars, however, would be 
coincident with star forming regions, therefore they can be excluded as the source of 
the majority of sGRBs. “Old” magnetar-like objects could also be formed through 
the accretion induced collapse (AIC) of an accreting WD or, maybe, as a result of a 
double white dwarf merger 1^221122^ For such old magnetars, however, no strong kicks 
are expected and therefore their distribution with respect to their host galaxy would 
differ from what is observed for sGRBs. It had also been suspectecG^HIHIIIl^ that 
nsns mergers may in some cases result in a rapidly spinning, high-held magnetar-like 
object and the recent limits on the maximum mass of neutron stars > 2 Me(2211211 
only make this a more intriguing possibility, see Sec. A magnetar engine would 
have a number of obvious advantages. There would be natural formation mecha¬ 
nisms and a long-lived or stable central object could plausibly be responsible for 
the observed late-time activity. A magnetar formed as a result of a nsns-merger 


’^Lee et find that nsbh collisions dominate by approximately a factor of five over nsns 

collisions. 
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would obviously inherit all the benefits of the nsns merger engine model. In ad¬ 
dition, having a magnetar at the engine of long and short bursts would naturally 
explain similarities between both types of bursts and with an energy emission rate 
(in the framework of the magnetic dipole model) of 

E^d oc (25) 

with B, R and oj being magnetic field strength, radius and angular frequency, such 
an engine could naturally produce the large variety of sGRBs that is observed. 

One could argue, though, that such a model, that depends on twelve powers of 
poorly known quantities has not much predictive power. Maybe more severe, there 
would be some challenge to explain both the prompt emission and the late-time 
activity, since the presence of a magnetar would initially produce an enormous 
baryon-loading which could prevent a burst to form in the first place. In fact, a 
recent studjl^ISl finds that collapse must occur before ~ 100 ms, otherwise a sGRB 
would be prevented. Maybe the choking of a GRB can be circumvented if the burst 
is only launched after a neutrino cooling time of several seconds when the neutrino- 
driven wind has ceased.l^^ 

Another alternative engine has been suggestecP^^ as a response to the discovery of 
late-time X-ray flares for short GRBs. In this model, the idea is that a neutron star 
accretes from a non-degenerate companion until it collapses into a bh surrounded 
by an accretion disk, generally considered the standard engine of GRBs. The X-ray 
flares would result when part of the relativistic ejecta interact with the extended 
companion star. In such a model, the neutron star would need to accrete ^ 0.7 Mq, 
not an entirely trivial task given that the Eddington accretion rate is MEdd ~ 10“® 
Mo/yr. Another question that would need further exploration is whether/for which 
EOS and rotation rate combination of the collapsing neutron star a substantial disk 
can form outside the ISCO of the forming black hole. 

In summary, there is likely room for several sGRB engines and maybe this would 
explain in part their diversity. However, while none of the suggestions is completely 
free of open questions, it is probably a fair statement that the best model to date 
for the bulk of sGRBs are compact binary mergers with -among them- a slight 
preference for nsns mergers, see the discussion in Sec. |5.3[ 

6. Summary 

The last years have seen a tremendous progress in our understanding of compact 
binary mergers, both on the theoretical and the observational side, and we have 
provided here an overview over various of their facets. 

The last decade has in particular witnessed the first detection of short GRB after¬ 
glows and subsequent observations have provided a wealth of information about the 
environments in which short GRBs occur. Most of the observed properties find a 
natural explanation in the compact binary merger model, but some properties such 
as very late activity keep posing a problem and are not understood. They could 
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either point to the merger event being way more complicated than imagined in to¬ 
day’s models or, alternatively, to engines that are different from the standard black 
hole plus disk picture. 

By now, practically all theoretical models agree that compact binary mergers eject 
enough material to be at least a major source of the heavy {A > 130) r-process and 
nucleosynthesis calculations show good agreement with the observed abundance dis¬ 
tributions in this regime. It has also become clear that a merger has several ways 
to enrich its host galaxy with neutron-rich matter: in addition to “dynamic ejecta” 
there are also neutrino- and/or magnetically driven winds and accretion tori that 
become unbound on viscous times scales and each of these channels has likely dif¬ 
ferent properties. There are even recent indications that the combination of the 
different channels may actually not just produce the heaviest elements, but even 
the whole r-process range, although contributions from other sources such as a su¬ 
pernovae are plausible, in particular for the lighter r-process elements. A question 
where no consensus exists yet, is whether compact binary mergers as dominant r- 
process source would be consistent with the chemical evolution of galaxies. While 
earlier work excluded them as dominant sources, more recent studies based on hy¬ 
drodynamic simulations come to the conclusion that compact binary mergers may 
well be consistent with the elemental scatter that is observed in stars of different 
ages. 

Maybe one of the most exiting new developments is the recent discovery of a 
credible radioactively powered, electromagnetic transient candidate, also known as 
’’macronova” or ’’kilonova”, that has been observed in the aftermath of a short GRB. 
In particular its time scale of about one week and the peak in the nIR are consistent 
with having been produced by very heavy (and therefore opaque) r-process mate¬ 
rial. If this is the correct interpretation, it connects for the first time directly short 
GRBs with compact binary mergers and r-process nucleosynthesis. Such transients 
will also increase the science returns in the era of gravitational wave astronomy that 
hopefully soon will begin. 
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